Universal intermittent properties of particle trajectories in highly turbulent flows 
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(Dated:) 

We present a collection of eight data sets, from state-of-the-art experiments and numerical sim- 
ulations on turbulent velocity statistics along particle trajectories obtained in different flows with 
Reynolds numbers in the range R\ £ [120 : 740]. Lagrangian structure functions from all data sets 
are found to collapse onto each other on a wide range of time lags, pointing towards the existence 
of a universal behaviour, within present statistical convergence, and calling for a unified theoretical 
description. Parisi-Frisch Multifractal theory, suitable extended to the dissipative scales and to the 
Lagrangian domain, is found to capture intermittency of velocity statistics over the whole three 
decades of temporal scales here investigated. 

PACS numbers: 



Understanding the statistical properties of particle 
tracers advected by turbulent flows is a challenging the- 
oretical and experimental problem [H, Q . It is a key in- 
gredient for the development of stochastic models [3, Q] , 
in such diverse contexts as turbulent combustion, indus- 
trial mixing, pollutant dispersion and cloud formation Q. 
The main difficulty of Lagrangian investigations, follow- 
ing particle trajectories, stems from the necessity to re- 
solve the wide range of time scales driving different parti- 
cle behaviours: from the longest, Tl, given by the stirring 
mechanism, to the shortest t v , typical of viscous dissipa- 
tion. Indeed the ratio, Tl/t ?) ~ R\, grows with the Tay- 
lor Reynolds number, R\, that varies up to few thousands 
in laboratory flows. Some aspects of Lagrangian statis- 
tics have been experimentally measured: particle acceler- 
ations , velocity fluctuations in the inertial range d, 0] 
and two-particle dispersion [1, Q. Others, connected to 
the entire range of motions, have long been restricted to 
numerical simulations [lj| EH, 0, 13, 14 1- A fundamen- 
tal open question is connected to intermittency, i.e. the 
observed strong deviations from Gaussian statistics, be- 



coming larger and larger when considering fluctuations at 
smaller and smaller scales. Besides, the dependency of 
velocity statistics at various temporal scales on large scale 
forcing and boundary conditions is the so-called problem 
of universality. Thus, universality features are linked to 
the degree of anisotropy and non-homogeneities of turbu- 
lent statistics (l5j. Similar problems have already been 



explored measuring the velocity fluctuations in the lab- 
oratory frame (Eulerian statistics), where clear evidence 
of universality have been obtained flil ] . 
To build a general theory of turbulent statistics, univer- 
sality is the first requirement and, if proved, ma y op en 
the possibility for effective stochastic modeling 17| in 
many applied situations. 

This Letter aims at investigating intermittency and uni- 
versality properties of velocity temporal fluctuations by 
quantitatively comparing data obtained from the most 
advanced laboratory @, @, Q and numerical 13, TI|, [H, 
Il3l . [Lsj experiments. Main outcomes of our analysis are 
twofold. First, we show that data collapse on a common 
functional form, providing evidence for universality of ve- 
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Rx 


Tr, (s) 


meas. vol. (v 3 ) 


Ntr 


Tech. 


Ref. 


1 


124 


8.5 x 10" 2 


340 3 


1.6 x 10 6 


PTV 


[8] 


2 


690 


9 x 10~ 4 


1700 3 


6.0 x 10 6 


PTV 


[7] 


3 


740 


2 x 10" 4 


6600 3 


9.5 x 10 3 


AD 


[6] 



TABLE I: Experiments. By columns: 1- experiment label; 
2- Taylor Reynolds number; 3- Kolmogorov time scale t v ; 4- 
measurement volume in unit of the Kolmogorov length scale 
r/; 5- N tr total number of Lagrangian trajectories measured; 
6- measurement technique: Particle Tracking Velocimetry 
(PTV) and Acoustic Doppler (AD); 7- Reference where infor- 
mation on the way the corresponding dataset was obtained 
can be found. 



locity fluctuations -up to moments currently achievable 
with high statistical accuracy. At intermediate and iner- 
tial scales, data show an intermittent behaviour. Second, 
we propose a stochastic phenomenological modelisation 
in the entire range of scales, using a Multifractal descrip- 
tion linking Eulerian and Lagrangian statistics. 

We analyse the probability distribution of velocity fluc- 
tuations at all scales, focusing on moments of these distri- 
butions, namely the Lagrangian Velocity Structure Func- 
tions (LVSF) of positive integer order p: 

S^(r) = ([v i (t + r)-v i (t)} p } = {(SrV i n (1) 

where i = x, y, z are the velocity components along a sin- 
gle particle path, and the average is defined over the en- 
semble of trajectories. As stationarity and homogeneity 
is assumed, moments of velocity increments only depend 
on the time lag r. In the inertial range, for r n «t«T;„ 
non-linear energy transfer governs the dynamics. Thus, 
from a dimensional viewpoint, only the scale r and the 
average energy dissipation rate for unit mass e should 
matter for the structure function behaviour. The only 
admissible choice is S^ p \t) ~ (er) p / 2 , but it does not 
take into account the fluctuating nature of energy dis- 
sipation. Empirical studies have indeed shown that the 
tails of the probability density functions of 5 T v become 
increasingly non-Gaussian at decreasing t/Tl. In terms 
of moments of the velocity fluctuations, intermittency 
reveals itself in the anomalous scaling exponents, i.e. a 
breakdown of the dimensional law for which we have that 
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with £(p) ^ p/2. Notice that when dissipative effects 
dominate, typically for scales r ~ r n and smaller, the 
power-law behaviour @ is no longer valid, and refined 
arguments have to be employed, as we will see in the 
following. 

The statistics of velocity fluctuations at varying time 
lag r can be quantitatively captured by the logarithmic 
derivatives of S^\t) versus sf \t) 0, 0, This 
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140 


256 3 


5 x 10 5 
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[11] 
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320 


1024 3 


5 x 10 6 
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[13] 
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400 


2048 3 


3 x 10 5 
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[10] 


4 


600 


1856 3 


1.6 x 10 7 
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[18] 
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650 


2048 3 


4 x 10 5 


N 


CS 


[12] 



TABLE II: Direct numerical simulations. By columns: 1- 
numerical simulation label; 2- Taylor Reynolds number Rx; 
3- number of collocation points iV 3 ; 4- total number of La- 
grangian tracers N p ; 5- characteristic of dissipation: normal 
viscous terms (N), weakly compressible code (C); 6- interpo- 
lation technique for Lagrangian integration: linear interpo- 
lation (L), tricubic interpolation (T); cubic-splines (CS); 7- 
Reference where information on the way the corresponding 
dataset was obtained can be found. 



defines the local scaling exponents 



d logger) 
dlogSfV) 
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For statistically isotropic turbulence, all components are 
equivalent, so that their spread quantifies the degree of 
anisotropy present in teh flow. The r-dependencc of 
(f>; T ) allows for a scale- by-scale characterisation of in- 
termittency 

Figure 1 shows the local exponents of order p ~ 4 from 
a collection of eight data sets, see Table I and II, for dif- 
ferent Reynolds numbers. Most of these data sets are 
new, as well as completely new is the performed anal- 
ysis, here presented. We focused on the fourth order 
moment, since it is the highest order achievable with sta- 
tistical convergence for all data sets. Two observations 
can be done. First, all data sets show a similar strong 
variation around the dissipative time t/t^ ~ 0(1) that 
depends on the Reynolds number, and then a clear ten- 
dency toward a plateau for larger lags r > 10t v . Second, 
all data sets, with comparable Reynolds numbers, well 
agree in the whole range of time lags. The relative scat- 
ter increases only for large t, due to the combined effects 
of the lack of statistics, the anisotropy of the flows and 
the different valus of R\. In particular, finite volume 
effects in experimental particle tracking can produce a 
small - but systematic - downward shift of the points at 
long-lag times [2l|, [22J. It is worth noticing that error 
bars estimated from anisotropic contributions decrease 
by going to small r, indicating that isotropy tends to 
be recovered at sufficiently small scales, i.e. large scale 
anisotropic contributions becomes less and less impor- 
tant. In addition, the fact that, at comparable Reynolds 
numbers, all data sets recover the same behaviour by 
going to smaller and smaller time lags provides a clear 
indication of Lagrangian universality of the energy cas- 
cade. Such an agreement has not been observed before 
and is comparable with that found for the corresponding 
Eulerian quantities (ljj . 
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FIG. 1: (Color online) Log-Lin plot of the fourth order local exponent, £(4,r), averaged over velocity components, as a 
function of the normalised time lag t/t v . Data sets come from three experiments (EXP) (see Table 1) and five direct numerical 
simulations (DNS) (see Table 2). Error bars are estimated from the spread between the three components, except in EXP3 where 
only two components were measured. Each data set is plotted only in the time range where recognised experimental/numerical 
limitations are not affecting the results. In particular, for each data set, the largest time lag always satisfies r < Tl. The 
minimal time lag is set by the highest fully resolved frequency. The shaded area displays the prediction obtained by the MF 
model by using DLo(h) or Drrih), with f3 = 4, for a range of R\ 6 [150 : 800], comparable with the range of R\ in the data. 
Notice that the MF predictions have been obtained by fixing equal to 7, the multiplicative constant in the definition of t v . 
The straight dashed line corresponds to the dimensional non-intermittent value ("(4, r) = 2, achieved at small time lags where 
structure functions do become differentiable. Notice that two among the DNS are sufficiently resolved to get the mentioned 
dimensional scaling in the high frequency limit. 



The quality of data shown in (Fig. 1) opens the pos- 
sibility to quantitatively test phcnomenological models 
for LVSF, scale-by-scale. Parisi-Frisch Multifractal (MF) 
model of the incrtial range 12311 . and its generalization 
to the dissipative range [24 



has proved to 



[251,126 

give a satisfactory description of Eulerian and Lagrangian 
fluctuations 



14L 28, \2 



fj. It is thus appealing to 
search for a link between Eulerian and Lagrangian statis- 
tics [13, 28, [29|, [3^1, since this points to a unique inter- 
pretation of turbulent fluctuations. Moreover, it would 
reduce the number of free parameters. According to 
the MF model, Eulerian velocity increments at inertial 
scales are characterised by a local Holder exponent h, i.e. 
5 r u ~ r h , whose probability V{h) ~ r 3 ~ D ^ is weighted 
by the Eulerian fractal dimension D{h) of the set where 
h is observed [23[. The dimensional relation r ~ r/S r u 
bridges Lagrangian fluctuations over a time la g r to the 
Eulerian ones at scale r. Following Refs. 27. |2&|. it is 
shown in Ref. [3(3] how to extend the MF framework to 
get a unified description at all time scales for Lagrangian 
turbulence. Accordingly, Lagrangian increments display 
a continuous and differentiable behaviour at the transi- 



tion from the dissipative to the inertial range, 

2h 





(4) 



[3 being a free parameter controlling the crossover around 
T~Tr), and Vq the root mean square velocity. In order to 
get a prediction for the behaviour of the LVSF, given by 



<(M P ) 



dhP h {T,T v )[5 T v{h)Y , 



(5) 



we have to consider, in ([4]), the intermittent fluctua- 
tions of the dissipative scale [3, 0, H^], T v (h)/TL ~ 
i?f" 1)/(m) . The last necessary ingredient is to spec- 
ify the probability of observing fluctuations of h. This is 
done in analogy to Eq. (g]): 



P h {T,T 11 )=Z- 1 {T) 



3-D(h) 



(0) 



where Z is a normalizing function [30( and D(h) the frac- 
tal dimension of the support of the exponents h. Once 
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specified the Reynolds number, we are left with two pa- 
rameters - the expoent (3 and a multiplicative constant 
in the definition of t v -, while the function D{h) comes 
from the knowledge of the Eulerian statistics. 
Eulerian Velocity Structure Functions (EVSF) have been 
measured in the last two decades (see Ref. [16[ for a data 
collection) providing a way to estimate the function D(h) 
based on empirical data. Many functional forms have 
been proposed in the literature [23| that are consistent 
with data, up to statistical uncertainties. Eulerian veloc- 
ity statistics can be measured in terms of longitudinal 
or transverse fluctuations. Fluid velocity along parti- 
cle paths is naturally sensitive to both kinds of fluctu- 
ations. We thus evaluated the LVSF in ((5]) using the 
fractal dimensions Dl {H) and Dxr{h) obtained by lon- 
gitudinal [lit an d transverse [3~0 ] moments of Eulerian 
luctuations, respectively. 

The shaded area in (Fig. 1) represents the range of 
variation of the MF prediction computed from Dl (K) 
or Dxrih), measured in the Eulerian statistics (see be- 
low), and at changing Reynolds numbers. This must be 
interpreted as our uncertainty. The prediction works very 
well: all data fall within the shaded area. The role of the 
parameters is clear. Changing (3 modifies the sharpness 
and shape of the dip region at t v - the larger (3 the more 
pronounced the dip; while, changing the multiplicative 
constant in the definition of t v has no effect on the curve 
shape, but it rigidly shifts the whole curve along the time 
axis. 

Increasing the Reynolds number R\, the flat region at 
large lags develops a longer plateau. In the limit R\ —> oo 
the MF model predicts C(4) ^ 1.71 from D Lo (h) and 
£(4) ~ 1.59 from DT r (h) statistics. 

For the Eulerian D(h), we used the following log- 
Poisson [23j functional form, 



D(h) 



3(h - h*) 
log(7) 



3(h*-h) \ 



3-d*. (7) 



Different couples of parameters, (h* , 7), have been chosen 
to fit longitudinal and transverse Eulerian fluctuations. 
The parameter d* = (1 — 3h*)/(l — 7) is fixed by imposing 
the exact relation for third order EVSF. For the longitu- 
dinal exponents 0], we used (h* Lo — 1/9, 7l = 2/3) 
[23^ . For the transverse exponents, we used (h^ r = 
1/9, 7Tr = 1/2) which fits the data in Ref. [3l[ (see [IH 
for details.) 

This comprehensive comparison of the best available 
experiments and direct numerical simulations provides 
strong evidence of the universality of Lagrangian statis- 
tics. One important open question is the effect of a 
mean flow, as in turbulent jets [33| and wall bounded 
turbulence, where strong persistence of anisotropy 
may break the recovery of small-scale universality. 
We showed that a Multifractal description is in good 
agreement with data, even in the dissipative range 
where intermittency is significantly increased. The 



Multifractal description captures the intermittency at 
all scales with only a few parameters, independent of 
the Reynolds number. This is the universal feature of 
Lagrangian turbulence revealed by this study. There 
exists a long debate on the statistical importance of 
vortex filaments around dissipative time and length 
scales 0, 34 1 . Simulations [lfl, HO, HI] show that the 
dip region for r ~ t v can be depleted/enhanced by 
decreasing/increasing the probability of particles being 
trapped in vortex filaments. The Multifractal model is 
able to capture the intermittency around t v with the 
help of the free parameter (3. Different values of j3 should 
then correspond to different statistical weights of vortex 
filaments along particle trajectories. 
Only further advances in both experimental techniques 
and numerical power will allow us to test the same 
questions here addressed also for the higher order 
statistics. 
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